Electromagnetic Absorption in an Anisotropic Layered Superconductor 
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We calculate the infrared conductivity tensor of a layered superconductor considering two different 
order parameter symmetries: strongly anisotropic s-wave with line nodes, and pure (d x 2_ y 2) d-wave. 
The calculations are performed within the quasiclassical theory of superconductivity and include the 
effects of non-magnetic scattering processes. We discuss to what extent measurements of the optical 
absorption can be relied upon to distinguish between these two order parameter symmetries. 
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Introduction. The symmetry of the superconduc- 
|ting order parameter in the high-T c cuprates is cur- 
rently the source of considerable scientific debate. A 
^number of experiments have suggested the presence of d- 

^"iwave pairing (see for example 1,2), while recent ARPES 

^^measurements 3 appear consistent with an anisotropic s- 

OQwave state. This is an important issue to resolve since 
the symmetry of the superconducting order parameter is 

QQa vital component in the understanding of the underlying 

£" — pairing mechanism at work in the high-T c systems. 

To date much theoretical effort has been directed to- 



^wards understanding the effects of strongly anisotropic 
Impairing on the thermodynamic and equilibrium proper- 
^\ties of a superconducting system. This work includes 
"^studies of the behavior of the order parameter and ex- 
C^citation spectrum in the presence of impurities 5-7 and 
^interfaces. 8,9 These considerations have revealed qualita- 
l tive differences between the properties of anisotropic s- 
^wave order parameters in comparison to rf-wave systems, 
^however there exists no unambiguous experimental veri- 
fication as yet. In this paper we seek to go beyond the 
consideration of equilibrium properties by calculating the 
frequency-dependent current response of an anisotropic 
layered superconductor to an externally applied electro- 
magnetic field. A study of the temperature dependent 
microwave conductivity has recently been put forth by 
Hirschfeld et al. 4 and Borkowski & Hirschfeld. 10 We de- 
monstrate that under certain circumstances the dynami- 
cal properties of anisotropic superconductors are highly 
sensitive to the detailed structure of the order parame- 
ter, and thus may serve as a reliable probe for the order 
parameter anisotropy. 

We perform our calculations within a microscopic mo- 
del for layered superconductors which incorporates the 
effects of non-magnetic scattering processes. Conside- 
ring two different order parameter symmetries (strongly 
anisotropic s-wave with lines of nodes, and pure (d x 2_ y 2) 
d-wave), we calculate both the in-plane and c-axis cur- 
rent response. We observe striking qualitative differences 
in the electrical current response of the different order pa- 
rameter symmetries which are strongly correlated to the 
degree of (non-magnetic) scattering in the system. These 



differences are not due just to the differences in the ex- 
citation spectrum, but rather arise from the formation 
of a band of optically active Andreev bound states at 
frequencies below the gap edge. 

Microscopic Model. We consider the microscopic 
model discussed in Refs. 11, 12 (the interlayer diffu- 
sion model), which is based on the quasiclassical theory 
of superconductivity. This model is characterized by 
an infinite periodic stack of incoherently coupled two- 
dimensional Fermi liquids. The in-plane transport is 
taken to be of the usual Fermi liquid type (i.e. me- 
diated by charged quasiparticles propagating coherently 
with an in-plane velocity vf). The interlayer transport, 
on the other hand, is diffusive in nature, originating 
from incoherent scattering processes. Interlayer scatte- 
ring may take place through several different types of 
scattering processes such as electron-electron, electron- 
phonon, electron-impurity, etc. This model should be 
appropriate for systems whose c-axis transport is of the 
SIS (superconductor-insulator-superconductor-. . . ) type. 
This type of behavior has recently been observed in a 
class of high-T c compounds, 13 suggesting that this mo- 
del may be appropriate for certain high-T c materials. 

A detailed description of the interlayer diffusion mo- 
del, along with the derivation of both the in-plane and 
c-axis frequency-dependent conductivity, has been given 
elsewhere 11,12 and will not be repeated here. Instead, we 
present only a brief summary of the key features of the 
model along with a discussion of the necessary phenome- 
nological parameters. 

For simplicity we assume an isotropic cylindrical Fermi 
surface which we parameterize by the angle <f>. Within 
the quasiclassical formulation of superconductivity, our 
model is completely defined by specifying the form of 
the scattering self-energy, S^. This quantity is conveni- 
ently written as the sum of an in-plane and two inter- 
plane parts, T,i = s| + + where I is the 
layer index. In-plane scattering is taken to be isotro- 
pic so that the scattering self-energy can be written as 
s| (e,t) = Citi(c,t), where ti is the single impurity t- 
matrix, 
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expressed in terms of the angular averaged single particle 
propagator, c/£. Here c; is the effective concentration of 
scattering centers, «o = «o 1 is an isotropic scattering po- 
tential, and Nf is the total density of states (per spin) at 
the Fermi energy. Following Buchholtz and Zwicknagl, 14 
we eliminate the parameters c; and uo in favor of an effec- 
tive normal-state scattering rate, l/Vy, and a normalized 
scattering cross-section, a. The normalized cross-section 
is a measure of the relative strength of the scattering and 
ranges from a = for weak scattering (Born limit), to 
(7=1 for resonant scattering (unitarity limit). 

We assume the c-axis coupling to be weak, and thus 
write the interlayer scattering self-energy in the Born ap- 
proximation, 
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The gauge operators, Ut,t±\, are defined in terms 
of an averaged interlayer vector potential A^ i±1 (t) 

by Ui,£±i(t) = exp[-ff ^-ti(i)^3], where is the 
layer spacing. The effective interlayer scattering life- 
time, Tj_(<f>,<f>'), is taken to be anisotropic. We des- 
cribe this anisotropy phenomenologically as Tj_(<f>,<f>') oc 
exp [—7 cos (<f> — <f>')]. The Fermi surface angles <f> and <f>' 
give the in-plane directions of the quasiparticle velocity 
before and after scattering to an adjacent layer. The pa- 
rameter 7 specifies to what degree the scattered electrons 
"remember" their initial momentum. Isotropic scattering 
corresponds to 7 = 0, while extreme forward scattering 
corresponds to 7 — ► 00. Since we neglect coherent trans- 
port along the c-axis (i.e. we set the Fermi velocity along 
the c-axis to zero), the interlayer scattering self-energy is 
the only source of interlayer coupling in the model. 

The interlayer diffusion model contains, as a minimal 
set, five material parameters: the transition tempera- 
ture, T c , the density of states at the Fermi level, Nf , the 
Fermi velocity, Vf , and the in-plane and inter-plane scat- 
tering lifetimes, ry and tj_. All of these quantities can 
be deduced from normal-state measurements. In order 
to accommodate anisotropic pairing, we have also intro- 
duced a normalized in-plane scattering cross-section, a, 
and an interlayer scattering anisotropy parameter, 7. 

The in-plane electrical current density is given in terms 
of the Keldysh component of the quasiparticle propaga- 



tor, (f t 
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by standard equations of Fermi liquid theory, 15 

de f d<f> 
Airi f 27r 



eN 



v f ((f))Tr{f 3 rjf((f);e,t)}, (3) 



where T3 is the third Pauli matrix. The microscopic ex- 
pression for the interlayer current density was derived for 
isotropic interlayer scattering in Ref. 11, and is generali- 
zed below to anisotropic scattering, 
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We compute the electrical conductivity by calculating 
and jl i + i in the presence of a weak electric field, and 
then reading off the appropriate coefficient. The proce- 
dure is rather involved and we thus refer the reader to 
Ref. 12 for the details of the calculation and the resulting 
expressions. 

Results. We discuss the electromagnetic absorption in 
the superconducting state for each of the following order 
parameter models: 



&Asw(<t>) 

A DW ((f>) 



A [l + cos(4^)]/2, 
A [l + 3cos (4</>)]/4, 
A cos (2(f>). 



(5) 
(6) 
(7) 



In our notation the subscripts ASW , XSW , and DW 
denote anisotropic s-wave, extended s-wave, and rf-wave, 
respectively. In terms of the irreducible representations 
of the (tetragonal) group, Aasw(4') an d Axsw(<t>) 
transform like the A\ g (identity) representation, while 
^dw(<I>) transforms like the B\ g (d x 2_ y 2) representation. 
All three order parameters possess nodes on the Fermi 
surface, but only Axsw(<f)) an d Aow(<f)) change sign. 

We calculate the supercurrents and associated conduc- 
tivities by numerically solving the quasiclassical trans- 
port equations. 11,12 These solutions must be carried out 
self-consistently for the order parameter amplitude, Ao, 
the scattering self-energy, X^, and the quasiclassical pro- 
pagator, gi. The conductivity can then be calculated 
from the analytic expressions derived in Ref. 12 by car- 
rying out the necessary integrations. 
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Figure 1. The in- plane infrared conductivity in the unitarity 
limit (<J = 1) for the ASW (dotted line), XSW (dashed line), 
and DW (solid line) models in units of the d.c. Drude conduc- 
tivity, erf, (0) = e 2 NfVf ry . The inset shows the corresponding 
excitation spectra. 



In figure 1 we show the real part of the in-plane 
frequency-dependent conductivity for the three order pa- 
rameter models presented above, together with the cor- 
responding excitation spectra (inset). These data are 
for a fairly clean system with a dimensionless scattering 
rate a = h/(2irT^k B T c o) = 0.02, where T c o is the transi- 
tion temperature in the absence of scattering. We con- 
sider here resonant scattering (a = 1), since it displays 
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the most striking features, and we take the temperature 
equal to zero. As was pointed out previously, 5,7 even a 
small amount of scattering opens a gap in the excitation 
spectrum of the ASW model, while the other two models 
display a finite density of states even at zero energy. The 
low energy enhancement in the density of states of the 
XSW and DW models can be interpreted as a band of 
optically active Andreev bound states. 6,12 These qualita- 
tive features are also represented in the absorption spec- 
trum. The ASW order parameter has a finite gap for 
any nonzero lifetime, ry, which results in a vanishing ab- 
sorption below a critical frequency u> < u> CT ~ l/ r ||- The 
XSW and DW order parameters still possess nodes, ho- 
wever, and exhibit a significantly increased absorption for 
hui< Ao. The enhanced absorption at low to has a width 
of the order of the cross-over energy e* ~ ^//iAo/4t|| , and 
comes from transitions within the bound band (resonant 
scattering), while the small feature at Tito ps Ao is a re- 
sult of transitions from the bound band to the gap edge. 
Figure 1 demonstrates that the low-frequency behavior 
of the conductivity is highly sensitive to sign changes in 
the order parameter, even though the density of states 
displays rather slight differences. 

Surprisingly, the ui^O absorption in the XSW model 
is actually larger than in the DW model. This difference 
can be quantitatively accounted for in terms of phase- 
space arguments. At T = one can show that u\ (u> — ► 0) is 
just proportional to (the total number of nodes in A(<^)) 
x (the slope of A((^) at the nodes) -1 . Hence, for our or- 
der parameter models, a^, w (u>—>0) ~ (4)x(2) _1 which 
is less than cjjf^yy (w->0) ~ (8)x(2 v / 2)" 1 by a factor of 
1 1 \/2 (plus terms of the order ~ %/ fry Ao] ) . This estimate 
is in quantitative agreement with the numerical result in 
figure 1. Our analysis can be set on a more rigorous foo- 
ting by noting that for rather clean systems (e* 2 <C A 2 ) 
at T = 0, the low frequency limit of the in-plane conduc- 
tivity can be written approximately as 



e* 2 cos 2 (^) 
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where A r (<^) is the real part of the scattering- 
renormalized order parameter at e = 0. We can represent 
A r ((f)) in a very general way by 



A r (<j>) = A [1 - (3 + 13 cos (2n<f>) + Ai] 



(9) 



where Ai ~ 1/ry is the real part of the off-diagonal con- 
tribution to the scattering self-energy. Note that all three 
models being considered here may be represented in this 
way. If A r ((f)) possesses nodes, then the largest contri- 
bution to the integral in equation (8) comes from the 
regions where A r (<^) ps 0. In this case we obtain the 
simple result: 
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which is valid as long as min[A r (^)] < — e*. It is inte- 
resting to note that our result does not depend on the 
value of the "symmetry parameter" n in equation (9). 



Equation (10) implies that the quantity Ao0i(w — *■ 0) 
is relatively independent of both the scattering lifetime, 
T||, and cross-section, a; these quantities only enter in- 
directly through the scattering self-energy piece, Ai. 
In fact, for the DW model (n = 1,(3 = 1), the off- 
diagonal scattering self-energy vanishes (Ai = 0) so that 
Aoc^w — ► 0) ~ e 2 Nf v 2 h / it , which is completely inde- 
pendent of the degree of scattering. 16 Figure 2 shows a 
series of absorption spectra for an XSW order parameter 
with a dimensionless scattering rate a = 0.1 for several 
different values of a. Note that the u> — ► limits for the 
conductivity all lie within ~ 10% of each other while the 
zero-energy values of the corresponding excitation spec- 
tra are very different. In the case of weak scattering, the 
limiting regime is only realized at very low frequencies, 
however this region attains an appreciable width for lar- 
ger cross-section values. It is important to point out that 
equation (10) quantitatively accounts for the ui^O limits 
in figure 2, and thus allows one, in principle, to obtain 
an estimate for the gap-anisotropy parameter (3 from a 
knowledge of the low-frequency absorption spectrum. 
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Figure 2. The in-plane conductivity in the XSW model, for 
u = 0.0 (dotted line), a = 0.5 (dashed line), and a = 1.0 
(solid line) in units of the d.c. Drude conductivity, f|j(0) = 
e 2 Nfvjr\y The lo — > limits are indicated by hollow circles. 
The inset shows the corresponding excitation spectra. 

A similar analysis of the low-frequency limit for the 
c-axis conductivity shows that this frequency range is 
dominated by the in-plane scattering properties. For bre- 
vity we restrict our analysis to the diffuse transmission 
limit (7 = 0) and the specular transmission limit (7 — ► 00) . 
Any finite 7- value will only change the absolute value of 
the result but not its physical behavior. For weakly cou- 
pled layers at zero-temperature, the c-axis conductivity 
has the general form: 



0) 



2e 2 rf 2 f d<f> f d<f>' A(<?i;0)A(<?i';0) 



2tt / 2tt 



t±(<M0 



(11) 



where N(<f>;0) is the angle resolved density of states at 
e = 0, and <f> and <f>' refer to Fermi surface positions in 
two adjacent planes. In the diffuse transmission limit, 
r£ 1 ((/), 4> r ) = r^ 1 and we obtain the simple result 



3 



of(w 0) 



NjTj_ \ /</> 



(12) 



where (• • denotes a Fermi surface average. This 
is what one expects for incoherent quasiparticle tunne- 
ling between two identical superconductors. Note that 
the zero-energy density of states, (N(<f>; 0))^, depends 
in a complicated way on the in-plane scattering para- 
meters T|| and a. In the specular transmission limit, 
r^ 1 (c/), <f)') = 2iTT^ 1 6((f) — <f)') and equation (11) becomes: 



^( w ^ )~ 2g2rf2 / 



(13) 



Following the same analysis as in the case of the in-plane 
conductivity, we find for specular transmission the result: 



(t^(lj 0) 
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which depends explicitly on the inter-plane scattering li- 
fetime, tj_, and implicitly on the in-plane parameters ry 
and a through e* and Ai. Equations (12) and (14) show 
that no universal behavior is expected in the inter-plane 
transport. In fact, the c-axis infrared absorption spec- 
trum is very nearly a direct map of the density of states 
in the layers (for Tito < Ao). Again one finds that the 
ASW model opens a gap for any finite in-plane lifetime 
(always assuming that tj_ ry). Note that in gene- 
ral the cross-over energy e* is quite different for different 
anisotropic pairing states as well as for weak and strong 
scattering, providing a way to distinguish between these 
various scenarios. 
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Figure 3. The c-axis infrared conductivity for our three or- 
der parameter models in units of the normal state value, 
<Tn =2e 2 N f d 2 /tl, for diffuse t ransmission. 

In figure 3 we show the c-axis conductivity in the dif- 
fuse transmission limit for the three different order para- 
meter models. For the XSW and DW models we have 
plotted results for both the Born and unitarity limits (the 
differences are very slight for the ASW model). Note 
that the XSW model no longer has a larger u> — ► li- 
mit than the DW case, nor do any of the models obey a 
universal low frequency limit. This is the behavior expec- 
ted from equation (12), since the zero-energy density of 



states in the XSW model is lower than that of the DW 
model for this choice of parameters. The situation is not 
significantly altered in the case of specular transmission, 
except for slight quantitative differences. We note, howe- 
ver, that the value of j has a much more profound effect 
on the coherent c-axis transport (i.e. Josephson tunne- 
ling). In fact, one finds that a finite j is necessary for 
the DW model to exhibit a finite Josephson effect. 17 

Conclusion. We have shown that anisotropic super- 
conductors with lines of nodes exhibit, at T —> and 
w^0, a strongly enhanced infrared absorption, and we 
have derived explicit expressions for its magnitude in 
terms of the order parameter anisotropy. The in-plane 
conductivity becomes universal for a pure ri-wave or- 
der parameter, 16 while an extended s-wave pairing state 
displays a nearly universal behavior in the clean limit. 
Conversely, the c-axis conductivity (for diffuse interlayer 
coupling) is not universal, but rather resembles the den- 
sity of states in the layers. Our results imply that one 
could, in principle, quantitatively ascertain the order pa- 
rameter anisotropy through either a study of the scaling 
of the ijj —> absorption with impurity concentration, or 
by searching for an impurity induced gap. 
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